High-frequency ultrasonic speckle velocimetry in sheared complex fluids 
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High-frequency ultrasonic pulses at 36 MHz are used to measure velocity profiles in a complex 
fluid sheared in the Couette geometry. Our technique is based on time-domain cross-correlation of 
ultrasonic speckle signals backscattered by the moving medium. Post-processing of acoustic data 
allows us to record a velocity profile in 0.02-2 s with a spatial resolution of 40 /tm over 1 mm. After 
a careful calibration using a Newtonian suspension, the technique is applied to a sheared lyotropic 
lamellar phase seeded with polystyrene spheres of diameter 3-10 /an. Time-averaged velocity profiles 
reveal the existence of inhomogeneous flows, with both wall slip and shear bands, in the vicinity 
of a shear-induced "layering" transition. Slow transient regimes and/or temporal fluctuations can 
also be resolved and exhibit complex spatio-temporal flow behaviors with sometimes more than two 
shear bands. 



PACS numbers: 43.58.+Z Acoustical measurements and instrumentation 
ing - 47.50.+d Non-Newtonian fluid flows - 83.60.Wc Flow instabilities 



43.60.+d Acoustic signal proccss- 



I. INTRODUCTION 

Accurate measurements of flow velocities are essential 
for exploring and understanding fluid dynamics. Such 
measurements raise important technical problems when- 
ever one needs them to be nonintrusive. Indeed, many 
applications where direct access to the fluid is impossi- 
ble perclude the use of local probes such as miniature 
piezoelectric pressure probes, hot films or hot wires. In 
other cases, the introduction of a probe inside the fluid 
may perturb the flow or even the fluid structure itself, 
and may lead to incorrect interpretations of the measu- 
rements. Thus, in general, much insight can be gained 
from some remote, noninvasive sensing of the flow field. 

In the following, we focus on the case of complex fluids 
under simple shear. A "complex" fluid may be charac- 
terized by the existence of a "mesoscopic" length located 
somewhere between the size of individual molecules and 
the size of the sample [J. For instance, this intermediate 
length scale is the particle diameter in a colloidal suspen- 
sion, the diameter of an oil droplet in an emulsion, or the 
radius of gyration of polymer coils. The existence of such 
a supramolecular organization (or "microstructure" ) can 
lead to very complicated behaviors under flow and to in- 
homogeneous velocity profiles under simple shear 0, • 
Since the microstructure of a complex fluid is very sen- 
sitive to local deformations, it is quite obvious that a 
nonintrusive technique is required in order to measure 
velocity profiles. 

The most popular nonintrusive techniques, namely 
Particle Imaging Velocimetry (PIV) or Laser Doppler Ve- 
locimetry (LDV), are based on the interaction between 
light and seeding particles following the flow 0, Q . How- 
ever, many complex fluids, such as emulsions, slurries, or 
pastes, may not be transparent enough to allow the use of 
PIV or LDV. Nuclear Magnetic Resonance (NMR) offers 
the possibility to image opaque media 0, ||| but requires 
the use of powerful magnets and remains expensive and 
tricky to set up. On the other hand, ultrasound appears 



as an efficient, cost-effective tool to measure velocity pro- 
files in a large range of fluids. 

In this paper, we adapt a speckle-based ultrasonic ve- 
locimetry technique 0, to complex fluid flows. Al- 
though based on the classical principle of backscattering 
by particles suspended in the flow, our technique brings 
an original contribution to both the fields of acoustical 
flow measurements and of complex fluid studies. Indeed, 
by using high-frequency pulses (frequencies larger than 
20 MHz), we show that we are able to measure veloc- 
ity profiles in complex fluids sheared between two plates 
separated by 1 mm with a spatial resolution of about 
40 /xm. Depending on the required accuracy, a full ve- 
locity profile can be obtained typically in 0.02 s to 2 s, 
which makes it possible to resolve transient regimes or 
temporal fluctuations of the flow. 

The paper is organized as follows. We first explain in 
more details why local velocity measurements are cru- 
cial to the understanding of complex fluids. Examples 
are given that show that inhomogeneous flows may oc- 
cur even in simple shear geometries and at small imposed 
stresses. The third section is devoted to a brief review 
of existing ultrasonic techniques for measuring flow ve- 
locities. We then present the electronic setup and the 
data analysis used for high-frequency ultrasonic speckle 
velocimetry (USV). Section deals with the calibration 
step necessary to obtain quantitative estimates of the ve- 
locity. In Section l*VTl high- frequency USV is applied to a 
particular complex fluid: an aqueous solution of surfac- 
tant. The ultrasonic data reveal the existence of inho- 
mogeneous flows, with both wall slip and shear bands, as 
well as complex spatio-temporal behaviors during tran- 
sient regimes. Finally, in light of the present results, the 
technique is compared to other nonintrusive tools used 
in the field of complex fluid flows. We emphasize on the 
promising spatio-temporal resolution of USV and discuss 
its possible future applications. 
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II. SHEAR FLOW OF COMPLEX FLUIDS: WHY 
ARE LOCAL MEASUREMENTS ESSENTIAL? 



A. Shearing a Newtonian fluid in the Couette 
geometry 
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where H is the height of the Couette cell. The viscosity r\ 
of the sample can thus be computed from the rheological 
data 7 and a by r\ = er/7, as long as the flow remains 
laminar and stationary. 




FIG. 1: (a) Sketch of the Couette cell used in our experiments, 
(b) Enlargement of the shear flow inside the gap between the 
two cylinders in the case of a Newtonian fluid below the onset 
of the Taylor-Couette instability. 



FIG. 2: Examples of orthoradial velocity profiles observed in 
a small gap Couette geometry, (a) Newtonian fluid sheared 
below the onset of the Taylor-Couette instability, (b) Wall 
slip in a complex fluid: apparent slippage results from the 
presence of thin lubricating layers near the walls, (c) Shear 
banding in a complex fluid near a shear-induced transition: 
two differently sheared bands coexist in the gap of the cell. 



It is well-known that shear flows of simple fluids are 
laminar at low velocities and may become unstable above 
a critical Reynolds number due to hydrodynamic insta- 
bilities HJ. A classical shear geometry is the Couette 
geometry in which the fluid is sheared between two con- 
centric cylinders as sketched on Fig. 2] Usually the inner 
cylinder of radius R x is rotating while the outer cylinder 
of radius R 2 remains fixed. The two cylinders are thus 
called "rotor" and "stator" respectively. The distance 
e = R 2 — Ri between the two walls is called the "gap" 
of the Couette cell. This Couette geometry will be used 
throughout the paper as well as the cylindrical coordi- 
nates (r,0,z). We will denote x = r — R 1 the distance 
from the inner cylinder and v = (v r ,vg,v z ) the velocity 
vector. 

A Couette device can work in two different modes: ei- 
ther a torque T is imposed on the rotor and the rotor 
angular velocity f2 is measured, or f2 is imposed and the 
torque T is measured. For a Newtonian fluid at small 
velocities, the flow is stationary and purely orthoradial 
v = (0, vg(r), 0). The tangential velocity v(x) = vg(r) 
decreases from the rotor velocity v a = R t Q, down to zero 
on the outer wall. In the small gap approximation i.e. 
when e <C Ri, the velocity profile is linear as shown in 
Fig. Ufa). The shear rate 7 (r#-component of the rate- 
of-strain tensor) and the shear stress a (tangential force 
per unit surface i.e. r#-component of the stress tensor) 
are then almost uniform across the gap and are simply 



B. Inhomogeneous flows in sheared complex fluids 

1. Apparent wall slip 

In complex fluids, the picture of the Couette flow may 
be radically different. For instance, when a complex 
fluid is confined between two plates and depending on 
the roughness properties of the plates, the fluid velocity 
close to the walls may strongly differ from that of the 
walls (see Fig. EJb)) : the fluid slips. Slippage is usually 
explained by the existence of two thin lubricating layers 
in which the fluid structure is very different from that 
in the bulk. For instance, in a colloidal suspension, elec- 
trostatic or stcric effects next to the walls may induce a 
depiction of the particles on a small distance, leaving a 
much less viscous fluid film at the walls and leading to 
this apparent wall slip jlfjj . In concentrated oil- in- water 
emulsions, wall slip is due to the presence of very thin 
water layers close to the walls. 

Even if the bulk fluid remains Newtonian, the veloc- 
ity profile between the two plates may be quite different 
from that expected in a uniformly sheared fluid that does 
not slip (see Fig.[2{a)). Wall slip in complex fluids raises 
important industrial problems and is often difficult to 
detect and to assess quantitatively. In particular, com- 
puting the viscosity from Eqs. and can be very 
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misleading since the effective shear rate in the bulk is 
sometimes orders of magnitude smaller than vq /e [Tol | . 

2. Hydrodynamic and elastic instabilities 

Moreover, even in the absence of wall slip and even for 
Newtonian fluids, above a critical rotation speed of the 
rotor, the laminar flow becomes unstable and vortices 
develop due to the centrifugal force (Taylor- Couette in- 
stability). Upon further increasing the rotation speed, 
this three-dimensional, inhomogeneous but still station- 
ary flow turns into a time-dependent turbulent flow 
The procedure based on Eqs. and J2J for measuring 7 
and a and computing the viscosity then completely fails 
because the flow is no longer laminar and stationary. 

About fifteen years ago, complex fluids such as entan- 
gled polymer solutions were shown to present inhomoge- 
neous flows even when sheared at rotor velocities far be- 
low the onset of the Taylor-Couette instability. These un- 
stable flows result from purely elastic instabilities |lllll2j| . 
Although due to the elasticity of the fluid rather than 
visco-inertial effects, the phenomenology of such instabil- 
ities is very similar to classical hydrodynamic instabilities 
but occur at very low Reynolds numbers. 



focused on the local characterization of the fluid struc- 
ture using birefringence neutron, x-ray or light scatter- 
ing and showed the coexistence of bands of different mi- 
crostructures in the vicinity of shear-induced transitions 

01311. 

As far as the flow field is concerned, rather few local 
velocity measurements are reported in the complex fluid 
literature The presence of inhomogeneous flows in 

some surfactant mixtures ("wormlike micelles") was first 
unveiled using NMR [25ll2r|. However, the existence of 
shear bands in the velocity profiles was firmly ascertained 
only recently using dynamic light scatterin g ( DLS) in 
heterodyne mode, a technique close to LDV |2fl,l28|| . 

From the various examples cited above, it is quite clear 
that local measurements are essential, for instance to pre- 
cisely follow the shear bands in both space and time. In 
general, a velocity profile with about 20 points accross 
a 1 mm gap, i.e. with a resolution of 50 /xm, will carry 
enough spatial information on the flow to allow quantita- 
tive conclusions on possible inhomogeneous flow profiles. 

III. MEASURING FLUID FLOWS WITH 
ULTRASOUND: A SHORT REVIEW 

A. Continuous wave systems 



3. Shear-induced structures and shear-banded flows 

Beside those inertial and elastic instabilities, more sub- 
tle effects may come into play that are due to the very 
nature of complex fluids. Indeed, application of shear to 
a complex fluid may completely change its microstruc- 
ture leading to a new shear-induced structure. Such a 
phenomenon may be called a shear-induced transition or 
a material instability in the sense that it results from 
a strong coupling between the shear flow and the mi- 
crostructure of the fluid [lflj . 

A material instability usually shows up on rheologi- 
cal data as a sudden jump in the viscosity when shear 
is increased. For instance, if the shear-induced struc- 
ture is much less viscous than the complex fluid at low 
shear, then a large drop of the visocity may be observed 
and the transition is called a shear-thinning transition. 
During such a transition, the fluid is assumed to phase- 
separate into two states: a shear-induced state that flows 
at a given shear rate and coexists with the old structure, 
thus giving rise to two shear bands E 13. The 
corresponding shear-banded velocity profile is sketched in 
Fig. He). 

In order to investigate those transitions, most previous 
studies have focused on global measurements such as the 
viscosity of the sample as a function of time [l?! ITsL IT^ | . 
Indeed, the time series rj(t) sometimes display complex 
temporal fluctuations on time scales characteristic of the 
microstructure rearrangements. However, such global 
measurements do not provide any information on the spa- 
tial structure of the flow H. Thus, other studies have 



Contrary to electromagnetic waves, the phase of acous- 
tic waves is easily accessible using piezoelectric transduc- 
ers. Moreover, ultrasound may propagate deeply into op- 
tically opaque media. These two properties have led to 
the spectacular development of ultrasonic imaging tech- 
niques since the early 1950s particularly in the biomedical 
domain. 

When an ultrasonic wave travels through biological tis- 
sues, it gets scattered by density and/or compressibility 
inhomogeneities [2*1 . If the scatterer moves in a flow, 
like for instance red blood cells inside blood vessels, its 
motion induces a Dopplcr shift in the wave frequency. 
Early systems for measuring blood velocities were based 
on the estimation of this Dopplcr frequency shift using a 
monochromatic ultrasonic wave Q . 

Unfortunately, even if these "true Dopplcr" systems 
helped to detect the occlusion of blood vessels, they did 
not offer any spatial resolution. Indeed, if one uses a 
continuous wave, one gets a very good resolution on the 
frequency shift (and thus on the velocity) but one also 
loses any temporal information on the echo arrival time 
(and thus on the scatterer position, as explained below). 
Continuous Doppler systems may thus be very efficient 
in the case of a single scatterer suspended in a flow whose 
position is measured independently |3*cj . 

B. Pulsed wave systems 

In order to discriminate between different scatterers in 
space and to measure velocity profiles, a solution is to 
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use short acoustic pulses. Indeed, if an ultrasonic pulse 
is sent by a transducer through a scattering medium, 
a series of backscattered echoes can be recorded on the 
same transducer. The precise nature of the backscattered 
(BS) signal will be discussed at length in Sec. lIVI At this 
point, it is only important to note that the arrival times 
tk of the echoes are directly linked to the position y^ of 
the scatterer along the acoustic beam by yk — c tk/2, 
where c is the effective sound speed and the factor 2 
accounts for the round trip from the transducer to the 
fc th inhomogeneity and back to the transducer. 

The idea is then to follow the scatterer motion through 
the evolution of the BS signals: between successive 
pulses, echoes will move along with the flow and their 
positions may be tracked in time. Although the Dop- 
pler effect is not used in pulsed wave systems, these are 
still mistcrmcd "pulsed Dopplcr systems" in the litera- 
ture and the technique is referred to as "ultrasonic Dop- 
pler velocimetry" (UDV) or "ultrasonic velocity profiler" 
(UVP) [""""J. This confusion may be due to the fact that, 
owing to the constraint of real-time display in commer- 
cial devices, BS signals arc sampled only at a few given 
depths and that these sampled signals arc analyzed in 
the frequency domain ||. 

UDV has been successfully applied to various classical 
problems in hydrodynamics such as the Taylor-Couette 
instability [32| . cylinder wakes |3^|, or magnetic fluid 
flows j34|. The fluid has to be seeded with small par- 
ticles that scatter ultrasound. Such acoustic "contrast 
agents" are assumed to follow the flow as Lagrangian 
tracers. The UVP has also been used fairly recently for 
in-line rhe olog ical studies of concentrated suspensions in 
pipe flows [35j] . 

If one drops the constraint of real-time display, the 
whole BS signals may be stored for post-processing. Di- 
rect tracking of the BS signals in the time domain is the 
subject of recent technical developments such as ultra- 
sonic speckle velocimetry which will be addressed in the 
next Section. 



applied to the measurement of tissue motion or blood 
flow 0, 111- They have also been adapted to the UVP 
system and tested on pipe flows to provide high time res- 
olution (of about 0.5 ms per profile) |u]]. Recently, one 
of us introduced a combination of a 2D echographic sys- 
tem and speckle tracking that we called 2D ultrasonic 
speckle velocimetry (USV) This system allowed us 
to image str ong vortical flows in Newtonian fluids in two 
dimensions |43j . 

Here we propose to adapt such a technique to complex 
fluids confined in a Couette cell of gap e ~ 1 mm. As 
explained in Section III Bl one needs to measure veloc- 
ity profiles with a spatial resolution of typically 50 /im. 
However, frequencies commonly used in commercial ul- 
trasonic velocimeters are in the range / = 1-10 MHz, 
and with a typical sound speed c ~ 1500 m.s -1 , this 
yields a resolution of about A = c n / f ~ 0.15-1.5 mm. 
Thus, conventional ultrasound does not permit very fine 
measurements and one has to turn to utrasound with 
frequencies larger than 20 MHz. Again, high frequencies 
have first attracted the biomedical community. High- 
resolution ultrasonic images have been recently obtained 
in dermatology, ophtalmology and stomatology 14*3, l45l o r 
for measuring blood flow in the microcirculation |46t l47j . 
To our knowledge, no such high-frequency measurements 
(above 20 MHz) have been reported so far in classical 
hydrodynamics nor in the field of complex fluids. 

Note that the resolution cannot be increased indefi- 
nitely by increasing the frequency because attenuation of 
sound waves sharply increases with frequency (as f 2 in 
water). Attenuation is due to the combination of thcrmo- 
viscous absorption and scattering. The signal-to-noise ra- 
tio rapidly deteriorates as / increases. We found that for 
our application to complex fluids, / = 36 MHz realizes 
a good compromise between resolution and attenuation 
over propagation distances of a couple of centimeters. Fi- 
nally, since transducer arrays are not available at such a 
high frequency, the technique will have to be restricted 
to one-dimensional measurements. 



IV. HIGH-FREQUENCY ULTRASONIC 
SPECKLE VELOCIMETRY 

A. Speckle tracking techniques 

When the moving medium contains a lot of scatterers 
per unit volume, the BS signal results from the inter- 
ferences of all the backscattered waves and appears as a 
complex, high-frequency signal called ultrasonic speckle 
(see Fig. [5] for a typical speckle signal and Sec. IIV Dl 
for more details, in particular about the possibility of 
multiple scattering). The development of fast digitiz- 
ers with large memories has allowed the full recording 
of backscattered BS signals for post-processing. Various 
time-domain algorithms for tracking the motion of ultra- 
sonic echoes have been proposed *3a, US HH, Hil, [""(j . So 
far, such speckle tracking techniques have mainly been 



B. High-frequency USV electronic system 

Figure |21 presents our USV electronic system. Focused 
ultrasonic pulses are generated by a PVDF piezo-polymcr 
immersion transducer of central frequency / = 36 MHz 
(Panametrics PI 50-2). The focal distance is 11.6 mm and 
the active element diameter is 6.3 mm. The transducer 
bandwidth is 11 MHz at -3 dB. The axial and lateral 
resolutions given by the manufacturer at -3 dB are 30 fim 
and 65 /im respectively and the depth of field is about 
1 mm. 

The transducer is controlled by a pulser-receiver 
unit (Panametrics 5900PR). The pulser generates 
220 V pulses with a rise time of about 1 ns. The pulse 
repetition frequency (PRF) is tunable from to 20 kHz. 
The receiver is equipped with a 200 MHz broadband 
amplifier of maximum voltage gain 54 dB as well as a 
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FIG. 3: Electronic system used for high-frequency USV. 



set of selectable high-pass and low-pass filters. BS sig- 
nals are sampled at f s = 500 MHz, stored on a high- 
speed PCI digitizer with 8 Mb on-board memory (Acqiris 
DP235), and later transferred to the host computer for 
post-processing. 



C. Experimental setup in the Couette geometry 

The shear flow is generated in a Plexiglas Couette cell 
with R t = 24 mm, R 2 = 25.07 mm, and H = 30 mm. 
The rotation of the inner cylinder is controlled by a stan- 
dard rheometer (TA Instruments AR1000). The whole 
cell is surrounded by water whose temperature is kept 
constant to within ±0.1°C. The thickness of the stator is 
2 mm everywhere except for a small rectangular window 
where the minimal thickness is 0.5 mm in order to avoid 
additional attenuation due to the Plexiglas. 




FIG. 4: Experimental setup for USV in complex fluids sheared 
in the Couette geometry. c , c p , and c w stand for the speed of 
sound in the complex fluid, in Plexiglas, and in water respec- 
tively. The inset shows a picture of the transducer and the 
stator as seen from above once the rotor has been removed. 



Ultrasonic pulses are incident on the stator with a 
given angle 9 ~ 20° relative to the normal to the window 



in the stator as sketched in Fig. 0] They travel through 
Plexiglas and enter the gap with an angle 9 that is given 
by the law of refraction. Since the sound speed of the 
working fluid may differ from that of water, and since 
the precise value of 9 depends on the exact arrangement 
of the acoustic beam relative the stator, 9 will be taken 
as an unknown until a careful calibration procedure is 
completed (see Section IV) | . 

Once inside the fluid, ultrasonic pulses get scattered 
by inhomogeneities that can be either naturally present 
(oil droplets in an emulsion for instance) or artificially 
introduced to enhance the acoustic contrast. The total 
round trip for a pulse travelling from the transducer to 
the rotor and back to the transducer lasts about 15 /is. 

The position of the transducer is tuned so that the 

1 mm gap lies into the focal spot in order to optimize 
the signal-to-noise ratio. Typical recorded BS signals are 
1000 point long, which corresponds to a transit time of 

2 /is or equivalently to a depth of roughly 1.5 mm inside 
the medium. Allowing for a fraction of a millimeter be- 
fore the stator and after the rotor helps to locate the two 
walls. Note that spurious reflections on water-Plcxiglas 
or Plexiglas-fluid interfaces are minimized by carefully 
choosing the angle 9 . In any case, even in the absence 
of interfaces, such an angle is necessary to get a non-zero 
projection of the velocity vector along the acoustic axis 
as explained below. 

D. Data analysis involved in USV 

1. Basic principle of USV 




FIG. 5: Speckle signal recorded in a 1 % wt. suspension of 
polystyrene spheres. The voltage V(t) from the transducer 
output is directly proportional to the pressure pit). The sig- 
nal was normalized by its maximum value V , which yields 
V/Vo = p/po where po is the maximum pressure. The sam- 
pling frequency is / s = 500 MHz and the recorded signal is 
1000 point long. The rotor velocity is v = 23.5 mm.s . The 
dotted rectangle denotes the time window used in Figs, ©and 
□ 

Figure [5] shows a typical signal backscattered by a di- 
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lute suspension of polystyrene spheres and normalized by 
the maximum amplitude. The noise level is about 10 mV 
in all our experiments and typical signal amplitudes arc 
in the range 50-200 mV corresponding to a signal-to- 
noise ratio of 10-30 dB. 

Under the assumption of single scattering, the speckle 
signal received between times t and t + At can be inter- 
preted as interferences coming from scatterers located be- 
tween y and y + Ay, where y = c t/2, Ay = c Q At/2, and 
y is the distance from the transducer along the acoustic 
beam. The approximation that no multiple scattering 
takes place is thus essential in order to use the echo- 
graphic conversion rule y = c a t/2. Multiple scattering 
is avoided by carefully controlling the amount and the 
properties of the scatterers (see Sec.lvTl. 

When the fluid is submitted to a shear flow, the in- 
terference pattern changes as the scatterers move along. 
Figure [B] shows two signals corresponding to two pulses 
separated by ST = 1 ms and plotted over 4 acoustic pe- 
riods. A global shift to the right of St ~ 3 ns is clearly 
visible. With c ~ 1500 m.s" 1 , this corresponds to 5y = 
c St/2 ~ 2 /^m and to a velocity v y = Sy/ST ~ 2 mm.s -1 
projected along the acoustic axis. 
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p T (t')p T+ s T (t' + r) , 



(3) 



where the fc th time window is centered around t^ = 
t s + kAt/2 and t s ~ 13.5 fis is a reference time that 
corresponds to the beginning of the gap. t s depends on 
experimental parameters such as the working tempera- 
ture and on the exact arrangement of the transducer rel- 
ative to the Couette cell (see Sec.|Vj|. The times t^ cor- 
respond to the centers of the various time windows over 
which the signals are split. They will be later converted 
into positions yk at which the velocity is measured using 
Vk = Co(tk —t s )/2. Such windows correspond to slices of 
width Ay = 2X ~ 80 fim along the acoustic axis sepa- 
rated by Ay/ 2 ~ 40 /im. 




FIG. 7: Correlation function Cfe(r) obtained by using Eq. © 
on the two signals of Fig.|S](«). The solid line corresponds to 
the parabolic interpolation Ck (t) and the dotted line indicates 
the position 8t ~ 2.7 ns of the maximum of Ck- 



FIG. 6: Zoom over 4 acoustic periods of the speckle signal 
p/po. Bullets (•) correspond to the speckle signal p T shown in 
Fig.0received after a pulse is sent at time T. Open circles (o) 
correspond to the signal p T +sr recorded over the same time 
window after a second pulse is sent 1 ms later at T + ST and 
when shear is applied. The motion of the scatterers result in 
a shift to the right of the signal i. e. the scatterers move away 
from the transducer. The rotor velocity is v = 23.5 mm.s -1 . 



2. Cross-correlation algorithm 

In order to get an accurate estimate of the time-shift St, 
we use the simple cross-correlation algorithm described 
below. Two consecutive signals p T andp T+ar recorded af- 
ter pulses sent at times T and T+ST are cross-correlated 
over small time windows of width At = 4£ ac ~ 0.11 fis, 
where i ac = 1/f is the acoustic period. More precisely, 
the following cross-correlation coefficient is computed ac- 



As shown in Fig. [7| the correlation function Cfe(r) 
is then interpolated around its main maximum by a 
parabola and the value St of r that maximizes this 
parabola is extracted: 



C k {5t) 



C k (r) 



(4) 



where Ck denotes the parabolic interpolation of Ck ■ Note 
that St depends both on the position tk of the time win- 
dow (fast "ultrasonic time" ) and on the time T for which 
the cross-correlation is performed (slow "pulse time"). 
Thus St will be noted as a function of the two variables 
tk and T below. 

Finally, it is straightforward to convert time-shifts St 
into velocities using: 



Co St(t k ,T) 
2 ST 



where 



Vk = c {t k - t s )/2 



(5) 



(G) 
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is the position of the center of the window along the 
acoustic axis. The origin y = is taken at the stator- 
fluid interface. Equation JHJ thus yields the projection 
v y of the velocity vector along the acoustic axis at the 
discrete positions yk and at time T. 

3. Limits of the technique 

The procedure based on the interpolation of Cfc(r) al- 
lows us to measure time-shifts as small as St ~ 0.2 ns 
that correspond to displacements Sy ~ 0.15 /im. Since 
the PRF / PRF = 1/ST can be made as small as desired, 
infinitesimal velocities could in principle be measured. 
However, at very small shear rates, decorrelation of the 
BS signals may occur due to Brownian motion of the scat- 
terers or to low-frequency mechanical vibrations in the 
experimental setup. Consequently, we will avoid using 
PRFs smaller than 10 Hz and we estimate the minimum 
measurable velocity at about 1 ^m.s" 1 . 

On the other hand, the maximum measurable velocity 
is reached when 5y ~ A ~ 40 /im. Indeed, displacements 
greater than the acoustic wavelength cannot be measured 
using the cross-correlation algorithm described above un- 
less phase unwrapping and continuity conditions are im- 
plemented, which will not be considered here j^. With 
the highest achievable PRF / PRF = 20 kHz, a 40 /im max- 
imum displacement means that velocities above 0.8 m.s -1 
may not be accessed using USV. Note however that this 
limiting range of [1 /zm.s -1 - 0.8 m.s -1 ] applies to the 
velocity v y projected along the acoustic axis. This range 
may be extended by varying the incidence angle 9 (see 
Sec. ESI and Eqs. ©-© below). 

Finally, Eq. (JSJ) shows that a velocity profile could in 
principle be obtained by cross-correlating only two suc- 
cessive BS signals i. e. within a time interval 5T that can 
be made very small depending of the PRF. In practice, 
as seen on Fig. [5] the speckle amplitude is never uni- 
form. Locally, destructive interferences or the absence of 
scatterers may lead to signal levels too small to be ana- 
lyzed (see the signal of Fig. [S] around t ~ 14.1 or 14.4 /is 
for instance). Time windows where the signal amplitude 
does not reach a given threshold (typically 20 % of the 
maximum amplitude) are thus left out of the analysis. 
To recover a full velocity profile, some averaging is then 
performed over several sucessive cross-correlations: 

Vy(yk) = {vy(y k ,T)) T . (7) 

Figure [S] shows twenty successive BS signals recorded 
in a sheared Newtonian suspension with / PRF = 1 kHz. 
The ultrasonic pulses leave the stator and enter the gap 
at t s ~ 13.5 /its and a small fixed echo corresponding to 
the rotor position may be seen at t r ~ 15.0 /its (as well 
as a stronger echo at about 15.2 [is that corresponds to 
a multiple reflection from the Plexiglas window in the 
stator). The slopes of the traces left by the echoes in 
this two-dimensional (t, T) diagram are inversely propor- 
tional to the local velocities according to Eq. (JSJ). Thus, 
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FIG. 8: BS signals corresponding to 20 pulses sent every 
ST = 1 ms. The normalized pressure p/po is coded in gray 
levels. In order to highlight the areas with strong speckle sig- 
nal, the colormap was saturated so that white (black resp.) 
corresponds to p/po > 0.4 (< —0.4 resp.). Time windows 
where max(|p/p j) < 0.2 are not taken into account in the 
data analysis (see text). The gap lies between t 3 ~ 13.5 /us 
and t r — 15.0 /is. The horizontal axis is the fast "ultrasonic 
time" t and the vertical axis is the slow "pulse time" T. The 
rotor velocity is v = 23.5 mm.s -1 . 

the signature of shear is rather clear: velocities increase 
from the stator, where the scatterers remain almost fixed, 
to the rotor. Moreover, the presence of moving echoes for 
times t > t r is not surprising: these echoes simply cor- 
respond to scattering of the wave reflected on the rotor. 
One can see that the whole range of times from t s to t r 
gets covered by some significant speckle signal in about 
20 ms in that case. 

If better accuracy is required, more averaging can be 
performed to improve the statistical convergence of the 
method (see Fig. ^5] below) . Depending on the trade- 
off between accuracy and temporal resolution, we found 
that USV allows us to obtain a full velocity profile every 
0.02-2 s. 



V. CALIBRATION PROCEDURE USING A 
NEWTONIAN FLUID 

A. Newtonian suspension of "contrast agents" 

In order to test and calibrate the experiment, much ef- 
fort was devoted to developing "contrast agents" that fea- 
ture good scattering properties for high-frequency USV. 
Indeed, too weak a scatterer yields a signal level too low 
to work with, whereas too strong a scatterer leads to mul- 
tiple scattering. When the scatterer is small compared to 
the acoustic wavelength (Raylcigh approximation), three 
parameters control the amount of scattering: the scat- 
terer diameter, its compressibility and its density [2flj . 
Moreover, if one wants the contrast agents to remain 
Lagrangian tracers of the flow, their density has to be 
matched to that of the fluid. 

We found that homemade polystyrene spheres of diam- 
eter a = 3-10 /jm were large and hard enough to scatter 
36 MHz pulses efficiently but soft enough to prevent mul- 
tiple scattering when diluted at a weight concentration 
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of 1 % in water. Such polystyrene spheres were obtained 
by polymerization following Refs. 0> C3- A 1 % wt. 
solution of benzoyl peroxide in divynylbenzene is first 
prepared. Then, a polydisperse emulsion is obtained by 
stirring 40 g of the previous mixture in 59 g of water and 
1 g of sodium dodecyl sulfate at room temperature. Fi- 
nally, the emulsion is polymerized at 90° C for two hours. 
The prepared polystyrene spheres were washed with wa- 
ter three times before drying. The 1 % wt. suspension of 
such polystyrene spheres is Newtonian and has the same 
viscosity as water. 

B. Measurement of c 

In order to use Eqs. JSJ and (JSJl, the sound speed c of 
the working fluid has to be measured. We used a classi- 
cal transmission setup consisting of two transducers fac- 
ing each other. One of them, the emitter, remains fixed 
while the other works as a receiver and can be moved 
by a computer-controlled actuator. Pulses are recorded 
for various displacements Sy of the receiver and averaged 
over 100 sweeps. The time-shifts St between the arrival 
times of the pulses are measured and the sound speed c 
is given by the slope of Sy vs. St (see Fig. 01. 




8t (jxs) 



FIG. 9: Measurement of the sound speed Co dbS the slope of 
Sy vs. St (see text) for the Newtonian 1 % wt. suspension of 
polystyrene spheres (•) and for the seeded lamellar phase of 
Sec. I VII (o). Both data sets were obtained at 32°C and Sy = 
corresponds to a separation of about 10 mm between the 
emitter and the receiver. Linear fits yield c = 1507±10 m.s -1 
for the Newtonian suspension (solid line) and c = 1495 ± 
10 m.s - for the lamellar phase (dotted line). In practice, 
those two fits are almost undistinguishable. Typical errorbars 
on St are about the size of the symbols. 

Due to technical constraints, the measurements of c 
were performed with 25 MHz transducers. Tempera- 



ture is controlled to within 0.1°C by a water circulation 
around the sample. Figure shows that c = 1500 m.s -1 
leads to a very good approximation of the data obtained 
in the suspension of polystyrene spheres at our working 
temperature of 32°C. This value of c is in good agree- 
ment with the measurements reported in the literature 
for the sound speed of water vs. temperature |29j . In the 
following, we will assume that this value does not vary 
significantly with frequency between 25 and 36 MHz and 
all the data will be taken at the same temperature of 
32°C. 

Finally, note that the measurement of c using the 
transmission of ultrasonic pulses between two transduc- 
ers also allows us to check that no significant multiple 
scattering takes place in our fluids by looking at the 
transmitted waveform : the received signal remains al- 
ways as short as the emitted pulse and no long-lasting 
echoes typical of multiply scattered waves were ever 
recorded during the measurements of c . 



C. Measurement of t s 

Another important parameter is t s that gives the po- 
sition of the stator-fluid interface and strongly depends 
on the experimental conditions. 
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FIG. 10: St I ST as a function of the position t of the cor- 
relation window obtained in the Newtonian suspension for 
t»o = 2.9 (•), 5.8 (o), 11.7 (■), 23.5 (a), and 47.0 mm.s" 1 (x). 
The solid lines are the best linear fits of the data for 13.55 < 
t < 14.95 [is. The dotted lines indicate the Plexiglas-fluid 
interfaces at the stator (left) and at the rotor (right). The 
results were averaged over 50 series of 20 pulses. Typical 
standard deviations are of the order of the marker size. Inset: 
blow-up of St I ST vs. t close to the stator. 

Figure ITUl presents the results of our speckle tracking 
algorithm averaged over 50 series of 20 pulses like that 
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shown in Fig. EI The rotor velocity v was varied between 
2.9 and 47.0 mm.s -1 . The linear behavior of the time- 
shift St with the position t of the correlation window is 
a direct signature of the uniform shear flow inside the 
Newtonian suspension. Indeed, all the velocities consid- 
ered here are far below the onset of the Taylor-Couette 
instability which occurs when v e 3 / 2 /vR^J 2 ~ 41 (where 
v ~ 10 -6 m 2 .s -1 is the kinematic viscosity of our cali- 
bration fluid) i.e. when v ~ 200 mm.s -1 . 

For our dilute suspension, no wall slip is expected and 
St should go to zero at the stator: 6t(t s ) = 0. By lin- 
early fitting St /ST vs. t for various values of v , different 
estimates of t s are obtained that yield an average value 
t a = 13.54 ± 0.02 fis. The uncertainty on t s is due to 
the presence of the stator which leads to a small fixed 
echo on the signals and to small edge effects around t s 
(see inset of Fig. HOfl . It corresponds to an uncertainty of 
about 15 /jm on the position of the stator-fluid interface. 



D. Conversion to orthoradial velocity profiles v(x) 

Inserting the values of c and t s found above in Eqs. (J3J) 
and jnj yields the velocity profile v y (y). Thus, so far, we 
obtained measurements of the velocity vector projected 
along the acoustic axis. In order to get some physically 
more relevant information, we will assume from now on 
that the radial velocity is zero (or at least negligible when 
compared to the orthoradial velocity). This is of course 
the case for our Newtonian suspension below the onset 
of the Taylor-Couette instability. In complex fluids, in- 
homogeneous flows with wall slip or shear banding are 
purely orthoradial (see Sec. Ill B|) so that the above ap- 
proximation remains valid as long as hydrodynamic or 
clastic instabilities do not occur. Note that in the geom- 
etry discussed here, USV measurements are not affected 
by a non-zero vertical component of the velocity field. 

Assuming that v = (0, v, 0) and using standard trigo- 
nometric relations (see Fig. 1110 . one gets: 



v(x) 



Ri 



R 2 sin 9 



v v(y) ■ 



where 



= vM 



2R 2 ycos6- R t . 



(8) 



(9) 



In the small gap approximation, Eqs. JSJ and © reduce 
to v(x) ~ Vy/ sin# and x ~ e — ysin#. 

Thus the last parameter to calibrate is the angle 9. 
This is done by looking for a value of 9 for which v(x) 
coincides with the velocity profile expected for a Newto- 
nian fluid: 



v(x) 



x 



R 2 
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FIG. 11: Conversion of v y (y) to v(x) in the Couette geometry. 
The flow is assumed to be purely orthoradial. x is the radial 
distance from the rotor whereas y is the distance from the 
stator-fluid interface along the acoustic axis. The gap width 
was exaggerated to reveal curvature effects accounted for by 
Eqs. © and ©. 
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FIG. 12: Velocity profiles v(x) obtained from Eqs. JHJ and © 
for the data of Fig. ESI with c = 1500 m.s -1 , t a = 13.54 /us, 
and 6 = 13.5°. Rotor velocities are v = 2.9 (•), 5.8 (o), 
11.7 (■), 23.5 (a), and 47.0 mm.s" 1 (x). The solid lines cor- 
respond to the Newtonian velocity profiles given by Eq. ESI 
with the above values of v . Inset: dimensionless data v/v 
vs. x/e. 



where the last term results from the small gap approxi- 
mation. 

Figure IT21 shows that 9 = 13.5° leads to very good re- 
sults for all the rotor velocities v considered here. Nor- 
malized data v/v„ vs. x/e show that the Newtonian ve- 
locity profiles u~u (l — x/e) are recovered for all our 
data with the set of parameters c , t s , and 9 inferred 
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from the above calibration (see inset of Fig. I12fl . This 
calibration method is very sensitive to small variations 
of 9 and allows its estimation to within 0.2°. Note that 9 
depends on both the temperature and the nature of the 
fluid under study through c . As long as the temperature 
remains the same, the new value 9' of the incidence angle 
in a fluid whose sound speed c' a differs from that of the 
calibration fluid is simply given by SnelPs law: 



^sin, 

Co 



(11) 



where c a and 8 are the parameters measured for the cal- 
ibration fluid. 



E. Statistical convergence of the USV 
measurements 

Finally, Fig. 1131 shows the value of the v y measured at 
two different positions along the acoustic beam as more 
and more averaging is performed. It is clear that con- 
vergence is reached when N ~ 300 pulses arc taken into 
account. The standard deviation of the measurements 
decreases from about 6 % for N ~ 100 to less than 3 % 
for N > 400. 
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FIG. 13: Convergence of the velocity estimates as a function 
of the total number of pulses N taken into account in the 
analysis. v y is measured at y — 0.6 (o) and 0.8 mm (•) for 
v = 23.5 mm.s -1 . 



VI. APPLICATION TO A COMPLEX FLUID 

A. The lyotropic lamellar phase under study 

Let us now turn to an application of high-frequency 
USV to a complex fluid. We chose to test USV in a 



lyotropic lamellar phase composed of SDS (6.5 % wt.), 
octanol (7.8 % wt.) and brine (20 g.L^ 1 NaCl). This sur- 
factant mixture is known to display a series of structural 
shear-induced transitions |50Ll5lll52l |. At equilibrium, the 
system is composed of surfactant bilayers stacked with a 
smectic distance of about 15 nm. Under shear flow, this 
lamellar phase was shown to form a closely packed as- 
sembly of multilamellar vesicles [HlJ. The diameter of 
these vesicles, commonly called "onions," is typically a 
few microns. 
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FIG. 14: Typical evolution of the viscosity r\ with the ro- 
tor velocity v a for the SDS/octanol/brine system. A shear- 
thinning transition occurs as Vo is increased (dashed line). 
The pictures show diffraction patterns observed in the low 
shear regime (left) and on the high shear branch (right): a 
laser beam crosses the cell along the r direction and a CCD 
camera collects scattered light in the plane (9, z) |5ljl . 



At low shear velocities, the diffraction pattern obser- 
ved with light scattering is an isotropic ring indicating 
the absence of any positional order. However, when the 
velocity v a is increased, Bragg peaks appear in the diffrac- 
tion pattern showing the apparition of long-range order 
in the vesicle positions. As sketched in Fig. EI this 
disordcr-to-ordcr transition corresponds to a sharp drop 
in the sample viscosity by a factor of about 3 jHJ ■ Such 
a shear-thinning transition is also observed in colloidal 
suspensions and is called a "layering" transition j5^. As 
explained in Sec. Ill Bl inhomogeneous velocity profiles 
are expected in the vicinity of this shear-induced transi- 
tion where the two states (ordered and disordered states) 
may coexist in the gap of a Couette cell. 

As for its ultrasonic behavior, the "onion" assembly 
does not scatter 36 MHz pulses. Thus, the lamellar phase 
was seeded with the polystyrene spheres used in Sec. [V] 
in order to get some speckle signal from the fluid. We 
checked that the "layering" transition is not significantly 
modified by the addition of such contrast agents. Fig- 
ure ED presents two pictures of the seeded lamellar phase 
as seen under crossed polarizers and after shearing: the 
granular texture characteristic of disordered multilamel- 
lar vesicles does not seem perturbed by the addition of 
polystyrene spheres at 1 % wt. 
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FIG. 15: Sheared lamellar phase seeded with polystyrene 
spheres of diameter 3-10 fim and observed under crossed po- 
larizers at two different magnifications: (a) x20 and (b) x50. 



B. Evidence for inhomogeneous flows 

As shown on Fig. |5J the sound speed of the seeded 
lamellar phase was measured to be c = 1495 ± 10 m.s -1 . 
Since the experiments and the calibration arc performed 
at the same temperature and since the sound speeds of 
both fluids are the same (within less than 1 %), we used 
c„ = 1500 m.s -1 , t s = 13.54 fjs and 8 = 13.5° to obtain 
the velocity profiles shown in Figs. 1161 and II 71 

The two main features of the velocity profiles are as 
follows: 

a. Wall slip. As clearly shown by the insets of 
Figs. El and El the velocity close to the walls is never 
equal to that of the walls. On the contrary, v(x = 0) may 
differ from v by more than 30 %. One can also notice 
that slip velocities are comparable at the rotor and at the 
stator for the smallest rotor velocity but become highly 
dissymmetric when v > 20 mm.s -1 . A detailed study of 
wall slip in this system is left for future work. 

b. Shear banding. Figure lltil also unveils the pres- 
ence of shear bands when v > 20 mm.s -1 . A highly 
sheared band nucleates close to the rotor and progres- 
sively fills the gap as v is increased (see Fig. El for 
the highest velocities). The shear rate in the high shear 
band is about three times larger than that in the low 
shear band. This ratio remains roughly constant as 
long as two bands coexist. For v ~ 70 mm.s -1 , the 
highly sheared region has invaded the whole gap. When 
v a > 70 mm.s -1 , the flow is laminar again and looks sim- 
ilar to that observed before the "layering" transition at 
v < 20 mm.s -1 . 

Simultaneously to the present ultrasonic measure- 
ments, DLS data obtained in our group have confirmed 
the existence of shear-banded flows and strong wall 
slip effects in the same lamellar phase without seeding 
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FIG. 16: Velocity profiles v(x) measured in the lamellar phase 
for rotor velocities v = 9.8 (•), 19.6 (o), 24.5 (■), and 
29.4 mm.s -1 (a). The data were time averaged over sev- 
eral profiles (3 to 7 depending on v a ) obtained from series of 
1000 pulses recorded in about 2 s. Error bars correspond to 
the standard deviation of these profiles and account mainly 
for temporal fluctuations of the flow. The inset shows the 
dimensionless data v/v vs. x/e. 



polystyrene particles [54| . This last observation shows a 
posteriori that our "contrast agents" do not perturb the 
flow and can be treated as simple Lagrangian tracers. 
Both USV and DLS data are in good qualitative agree- 
ment with classical theoretical pictures of shear-banded 
flows 0, |5f| . Our purpose here is only to show the 
relevance of high-frequency USV in the field of sheared 
complex fluids and a more quantitative analysis of shear 
banding in light of DLS results can be found in Ref. [54| • 



C. Dynamical measurements during a transient 

Finally let us show how USV may be used to follow 
the dynamics of a complex fluid under shear. Until now, 
we have considered time-averaged velocity profiles. How- 
ever, during the "layering" transition, noticeable fluctua- 
tions of the flow field occur and result in large error bars 
for v Q = 24.5 or 29.4 mm.s -1 for instance. Actually, some 
very rich dynamics were reported in this sheared lamel- 
lar phase when the torque applied on the rotor is imposed 
instead of the velocity v [l9j. In that case, the viscos- 
ity r/(t) displays large, complex fluctuations on long time 
scales of the order of 100 s. Similar (although faster) 
temporal behaviors were also observed in other surfac- 
tant systems 0, 0] . 

The origin of these fluctuations and their precise sta- 
tistical properties remain unclear. One possibility is that 
the flow presents shear bands that may couple to each 
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FIG. 17: Same as FigED but with 
68.5 (■), and 97.9 mm.s -1 (a). 



v = 34.3 (•), 53.8 (o), 



other and lead to the observed behavior of the viscosity. 
Thus the dynamics could be not only temporal as obser- 
ved on T)(t) but also spatial as suggested in Ref. ppf. 

Figure ITS! shows a transient regime recorded after the 
torque applied to the rotor was suddenly increased from 
83 to 87 mN at time t = 0. The response of the rotor 
velocity v (t) is shown in Fig. I18f a~) [5g. A maximum 
of about 38 mm.s -1 is reached at t ~ 300 s before v (t) 
slowly relaxes with small fluctuations. Velocity profiles 
recorded during that transient show that: (i) during the 
initial increase of v (t) , a shear-banded flow is generated 
with three bands: two highly sheared bands at both the 
rotor and the stator and a low shear region in between 
(see Fig. llHf b)-(d)): (ii) the band located at the stator 
disappears when V (t) decreases, leaving only two bands 
similar to Figs.lrplandnTlfsee Fig.Oe)-(f)); (hi) states 
where the highly sheared band almost disappears can be 
observed (see Fig. EUg)) even though the rotor velocity 
v > 28 mm.s -1 indicates that a band should extend over 
almost half the gap if the average picture at imposed v 
remained valid when the torque is imposed (see the a 
symbols corresponding to v a = 29 A mm.s -1 in Fig. I1GII . 



VII. DISCUSSION AND CONCLUSIONS 
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FIG. 18: Transient regime observed when the torque is in- 
creased from 83 to 87 mN at t = 0. (a) Velocity v of the 
rotor as a function of time t. Bullets (•) correspond to the 
points where the velocity profiles shown in (b)-(g) were mea- 
sured: (b) t = 45 s, (c) t = 60 s, (d) t = 107 s, (e) t = 267 s, 
(f) t = 647 s, and (g) t = 1036 s. The velocity profiles were 
obtained from series of 800 pulses recorded in 1.6 s. 



A. Summary 



The above results show that high-frequency USV is 
a powerful tool for measuring velocity profiles in sheared 
complex fluids. We have successfully tested the technique 
on a lyotropic lamellar phase that undergoes a shear- 
induced "layering" transition. 

When the rotor velocity, i.e. the global shear rate, is 
imposed, shear banding is observed: time-averaged veloc- 



ity profiles display a highly sheared region that nucleates 
from the rotor and progressively invades the gap as v Q 
is increased. Wall slip can also be very large especially 
during the transition, which makes global measurements 
such as the sample viscosity tricky to analyze. 

Moreover, USV allowed us to record velocity profiles 
in 1 s typically. This acquisition time is short enough 
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(compared to the intrinsic time scales of our complex 
fluid) to follow the velocity profiles in time. We showed 
the existence of more complex, three-band flows during 
a transient at imposed torque i.e. imposed shear stress. 
This raises the question of the structure of the flow dur- 
ing the spontaneous oscillations of the viscosity observed 
in the same lamellar phase 0, [2(| • USV may thus pro- 
vide a "time-resolved" tool for investigating spatially the 
dynamical behaviors of r](t). 

B. Comparison with other local techniques 

In this paper, we have shown that we were able to mea- 
sure velocity profiles with a spatial resolution of about 
40 yum and every 0.02-2 s. High-frequency USV com- 
pares very well with other local techniques for measuring 
flow fields. 

LDV, or cquivalently DLS, provides in principle a sli- 
ghtly better spatial resolution of 10-50 fiva but only 
yields a pointlikc measurement. In order to obtain a full 
velocity profile with LDV or DLS, one has to mechan- 
ically scan the whole gap of the cell [27j]. This usually 
takes about 1 min and may raise major interpretation 
problems when the flow evolves on time scales shorter 
than 1 min. 

NMR appears as a very promising technique for com- 
plex fluids |2f| |2(| since no seeding is required and 
two-dimensional images of the flow are obtained. How- 
ever, acquisition times are rather long too (from about 
20 s to more than an hour depending on the number of 
points in the image), making it rather difficult to access 
dynamical local information. 

The main drawback of USV is that it needs a good con- 
trol of the scattering properties of the fluid and that it 
may require that the flow be seeded with contrast agents. 
Unless their mechanical properties is strictly controlled, 
most concentrated suspensions or emulsions will lead to 



multiple scattering of high-frequency pulses. Thus, fur- 
ther application of USV to other complex fluids will in- 
volve a careful study of their acoustic properties. How- 
ever, we believe that USV is likely to extend the range 
of complex fluids in which local measurements may be 
performed, in particular to optically opaque fluids. 



C. Perspectives 

The biggest advantage of USV lies undoubtedly in its 
ability to perform nonintrusive velocity measurements 
fast enough to follow the dynamics of complex fluids un- 
der low shear. In principle, in an ideal scattering envi- 
ronment, the temporal resolution of USV could even get 
as short as 1 ms per measurement. 

Future research directions on USV are as follows: (i) 
more measurements on seeded lamellar phases will be 
performed and compared to DLS data; more precisely, 
we will focus on the spatio-temporal flow dynamics in 
regimes where sustained oscillations of the viscosity are 
observed; (ii) we will investigate shear flows of other flu- 
ids such as emulsions, gels, or granular pastes; (iii) we 
also intend to apply high-frequency USV to the case of 
very small deformations in order to perform some "lo- 
cal linear rhcology" and understand how a complex fluid 
starts to flow in the transition from linear to nonlinear 
rheology. 
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